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Abstract 

We present a list of post-AGB stars showing molecular line absorption 
and emission in the optical spectrum. Two objects show CH + , one in emis- 
sion and one in absorption, and ten stars show C2 and CN in absorption. 
The Doppler velocities of the C2 lines and the rotational temperatures 
indicate that the line forming region is the AGB remnant. An analysis 
of the post-AGB stars of which CO millimeter data is available suggests 
that the C2 expansion velocity is of the same order as the CO expansion 
velocity. HD 56126 has been studied in detail and we find a mass-loss 
rate of M = 2.8 x 10" 4 Meyr -1 , /c 2 = 2.4 x 1(T 8 and / CN = 1.3 x 10" 8 . 
The mass loss derived from C2 is significantly larger than M = 1.2 x 10 -5 
M©yr _1 derived from CO. We find that all objects with the 21/im fea- 
ture in emission show C2 and CN absorption, but not all objects with C2 
and CN detections show a 21/^m feature. 

keywords: molecules - physical conditions in AGB remnant - mass-loss 
history on AGB 

'Based on observations with the WHT/UES 
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1 Introduction 



The study of molecular lines in the optical spectra of post-AGB stars started 
with the paper by Waelkens et at (1992) in which they presented a band 
of narrow emission lines in the spectrum of the Red Rectangle (HD 44179). 
These emission lines were identified as the (0,0) C3 band of CH + by Balm 
and Jura (1992). The presence of the Swan and Phillips bands of C2 (Fig. 
and red system bands of CN (Fig. [2]) in the spectrum of HD 56126 studied by 
Bakker (1994), showed that the outflow velocity and the excitation conditions 
of C2 suggest that these lines are formed in the AGB remnant. This opens 
new possibilities to study the physical conditions in AGB remnants. Hrivnak 
(1995) presents a list of post-AGB stars showing C2 and C3 in absorption. Here 
we present the observations of post-AGB stars that show molecular lines in 
the optical spectrum. For several stars the excitation temperature and outflow 
velocity have been determined. The C2 absorption in HD 56126 is studied in 
detail by modeling the population density over the rotational quantum number 
(Section |). 

2 Detection of Molecules 

In Table we list the first results of our work. For each star the 3th and 4th 
column give the expansion velocity derived from the optical (C2,CN,CH + ) and 
the CO millimeter or OH maser line respectively. The 5th and 6th column 
give the rotational temperature and column density derived from C2 (3,0) or 
CH + (0,0) band. Negative velocities are due to outflow of material. The C2 
line absorption is formed closer to the star than the CO line emission, the two 
different molecules trace different material and hence different stages of the 
AGB evolution. Our data may indicate that the expansion velocity decreases 
as the star evolves along the AGB (Table [j]) . However a larger sample is needed 
to quantify this result. Stars showing the unidentified 21/im emission feature 
all exhibit C2 and CN absorption, strongly suggesting that the 21/nm feature 
is from carbon-rich material, but not all stars with C2 absorption show the 
21/im feature. 

3 Modeling the C2 absorption Bands of HD 56126 

From the rotational diagram of the C2 bands (Fig. [|) we find T rot = 240 K and 
for the CN (1,0) band we find T mt = 24 K with a column density of Nc 2 = 
20 x 10 14 cm" 2 and A CN = 11 X 10 14 cm" 2 . Unlike CN, C 2 is a homo-nuclear 
molecule having no permanent dipole moment and the rotation temperature 
is not a good indicator of the kinetic temperature because of pumping by the 
stellar radiation field. We use the CN rotational temperature as a measure of 
the kinetic temperature of the gas. 
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Tabic 1 : Expansion velocities of the AGB remnant 



Object 


Id* 


^optical 


Sv+ 


T 

rmrot 


A^Boltz x 10 14 


Remark 








[km s~ 




[K] 


[cm- 2 ] 






IRAS 04296+3429 


c 2 


-7.5 


-12 


138 


14 






IRAS 05113+1347 


c 2 


-5.2 


no 


198 


12 


not observed in 


CO 


HD 56126 


c 2 


-8.8 


-10.0 


240 


20 






IRAS 08005-2356 


c 2 


38.1 


-50.0 


150 


31 


OH maser 




IRC +10216 


c 2 


17.9 


-15.7 






AGB star 




AFGL 2688 


c 2 


15.0 


-22.8 


56 


26 






HD 44179 


CH+ 


—4.5 em. 


-3 


267 


2 






HD 213985 


CH+ 


-5.4 abs. 


nd 






not detected in 


CO 



* all stars with C 2 do also show CN absorption 

+ from CO millimeter emission. For IRAS 08005-2356 there is no CO data 
available and we used the velocities from the OH maser line 
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Figure 1: The Phillips (2,0) band in the spectrum of HD 56126. The C 2 ab- 
sorption lines are not resolved. 
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Figure 2: The CN A 2 II - X 2 E+ (1,0) band in the spectrum of HD 56126. The 
computed synthetic spectrum is shifted down by 0.5. 
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Rotational Diagram of the C 2 A 1 ^ - X 1 E+ 

Figure 3: (1,0), (2,0), and (3,0) absorption band. The P, Q, and R branches 
are denoted by a plus, asterisk, and a square respectively. Only the weakest 
rotational band (3,0) is optically thin and gives an almost linear relation in the 
rotational diagram with T rot = 240 K and N Bo \tz = 20 x 10 14 cm~ 2 . 
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Figure 4: Relative rotational diagram of the C 2 A 1 !^ - X 1 S+ (1,0), (2,0), 
and (3,0) absorption band of HD 56126 with the model fit for T^ n = 24 K 
superimposed. 
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Table 2: Overview of derived mass-loss rates for HD 56126 
Modeling C2 excitation 2.8 x 1CT 4 M Q yr _1 this study 

Infrared excess from dust 1.9 x 10~ 4 M Q yr _1 this study; gas/dust=100 
CO millimeter emission 1.2 x 1CT 5 Meyr^ 1 12 CO(2-l) Omont et al. 1993 

We have tried to fit the relative rotational diagram (Fig. ^) with the pop- 
ulation density distribution derived from modeling the C2 excitation by taking 
into account radiative pumping and collisional de-excitation (Van Dishocck and 
Black 1982). The best fit is reached for n c a/I = 3.25 x 1CT 14 . Taking the latest 
H2-C2 cross section of a = 7.8 x 1CP 16 cm -2 (Phillips 1994) the mass-loss rate 
can be calculated using Eq. [IJ This gives a mass-loss rate of Mc 2 = 2.8 x 10 
Meyr" 1 . 

„> oc , n - 4 ( v ex P n c a 4.1xl0 5 \ _ x 

Fitting an optically thin dust model (Waters et al. 1988) to the the spectral 
energy distribution of HD 56126 yields a dust inner radius of 2 x 10 3 i?* and 
M = 1.9x 10~ 4 M©yr _1 . Assuming that the molecular line absorption originates 
at this dust inner radius we find a particle abundance relative to H2 of fc 2 = 
2.4 x 10~ 8 and /cn = 1.3xl0~ 8 . Table ^ summarizes the mass- loss rates derived 
for HD 56126, where we have scaled the different mass-loss rates to i?* = 50i?o, 
T e s — 6500 K (logL = 3.6) and D = 2.7 kpc. The mass-loss rate derived from 
C2 and from the IR excess are significantly higher than that derived from the 
CO millimeter emission, which might indicate that the mass-loss rate increased 
dramatically towards the end of the AGB. 
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